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I. INTRODUCTION 



The wide application of functional determinants in Quantum and Statistical Physics is by now a well known fact. 

■ Typically, one is faced to the necessity of defining a regularized determinant for elliptic differential operators, among 
[ which the Dirac first order one plays a central role. An interesting related problem is the modification of physical 

04 ■ quantities due to de presence of boundaries. The study of boundary effects has lately received much attention [Q-^, 
^ ' since it is of importance in many different situations, like effective models for strong interactions , quantum cosmology 
and application of OFT to statistical systems , among others. 

In previous work [ p^Jll| , we have studied elliptic Dirac boundary problems in the case of local boundary conditions. 
, In particular, we have developed for this case a scheme for evaluating determinants from the knowledge of the 

■ associated Green's function, based on Seeley's theory of complex powers . 

\^ ' Another type of boundary conditions extensively studied in the literature are global ones, of the type introduced 



On by Atiyah, Patodi and Singer (APS) 13 1 in connection with the index theorem for manifolds with boundaries (see 
[ p^ for a review.) Other motivations for considering these global (or spectral) conditions are their consistency with 
charge conjugation and chiral invariance, and the presence of topological obstructions for the chiral Dirac operator 
O ^ under local boundary conditions (although this restriction no longer holds when considering the whole Dirac operator 

0^. 0.) , 

i-G ' In this paper we present the complete evaluation of the determinant of the Dirac operator on a disk, in the presence 
of an axially symmetric Abelian flux, under spectral boundary conditions (see [ p5| , p^ and references therein for related 
•'-H . work). 

' In section II we establish our conventions and set the problem to be considered. 
^ . In Section III we study zero modes and obtain the projector on the null space of the operator. 

Section IV is devoted to the evaluation of the Green's function of the problem in the subspace orthogonal to such 
null space. The knowledge of this Green's function is necessary in the evaluation of the determinant, which is done 
in Section V. This evaluation is performed in two steps. The first one involves the consideration of the quotient of 
determinants with and without gauge field, under a fixed boundary condition. For this situation, we grow the field 
along a continuous path through a parameter a, and perform a point splitting regularization of the a-derivative of 
the logarithm of the quotient of determinants, written in terms of the Green's function. 

The second step amounts to the evaluation, via ^-function, of the quotient of the free operators, under different global 
boundary conditions, which is possible thanks to the explicit determination of their eigenfunctions and corresponding 
eigenvalues. 

In Section VI we discuss the application of the APS index theorem to our example, with emphasis on the effect of 
the boundary isomorphism a appearing in the Dirac operator. 
Finally, Section VII contains our conclusions. 



*Partially supported by CONICET, Argentina. 



1 



II. SETTING OF THE PROBLEM 



We will evaluate the determinant of the operator D =i^+ 4 acting on functions defined on a two dimensional disk 
of radius i?, under APS boundary conditions. 

We take to be an Abelian gauge field. As it is well known, it can be written as = e^^ + S^jy (eoi = 
— eio = 1). We set ?7 = 0, thus choosing the Lorentz gauge. For cf) we take a smooth bounded function (j) = 4>{^)\ then 
Ar = Q and Ae{r) = -^^^(r) = -cj) {r). We call 

K=^^^<j> Ae R de = -R4>' (R). (1) 
Our convention for two dimensional Dirac matrices is 

I) ' ^i = (z V) ' ^'^{l -l)' (2) 

which satisfy 

l/ili^ =Sn„ I +i tfiv 75- (3) 
Therefore, the free Dirac operator can be written as 

i ^ = i (^0 ao +71 = 2z ^ ^ , (4) 

where X = x^^-Vi x\ and X* = xq — i xi or, in polar coordinates 



with 



i ^ = ii-jr dr + -7e dg), (5) 



With these conventions, the full Dirac operator can be written as 

D = e-^^C-) i ^ e-^C-) = ( ^ ^ '^''^^ + ^) ) , (7) 

with 

a = - i (8) 

and 

B{r) = -- de+ dr<P{r)- (9) 
r 

At the boundary, the self adjoint operator B{r) becomes 

B = B{R) = ~[idg + Kl (10) 

with eigenvectors e'"^ and eigenvalues -j^(ri — k). We take the radial variable to be conveniently adimensionalized 
throught multiplication by a fixed constant with dimensions of mass. 

We will consider the action of the differential operator D on the space of functions satisfying global boundary 
conditions, characterized by 

with 



2 



V> = ^ .), (12) 



2tt 

n>k+l 



where k is the integer such that k < K < k + 1. Notice that 

a{l-V>)a* ^aV^a* ^ e™« (e-^ • ) = 7'< (13) 

n<k+l 

and the operator so defined, {D)^, turns out to be self adjoint. The presence of a as in (^) has a relevant consequence 
on the boundary conditions for the lower components. As we will discuss later, this also reflects on the form of the 
index theorem for the present situation. 

Our aim is to compute the quotient of the determinants of the operators {D)i^ and {i ^)k=o- Since the global 
boundary conditions in Eq. (^) depend on the flux $ as a step function, there is no continuous family connecting 
both operators. So, we will proceed in two steps: 

{D)n - -> ^)Ai=o. (14) 

In the first step, where there is no change of boundary conditions, we can grow the gauge field by varying a from 
to 1 in 

= i ^ + a / = 6""'^=^'^^ i ^ e-a750(r)^ ^-^^-^ < a < 1, (15) 

thus going smoothly from the free to the full Dirac operator. The explicit knowledge of the Green's function will allow 
us to perform the calculation of this step, where we will use a gauge invariant point splitting regularization of the 
a-derivative of the determinant. The second step will be achieved by using a ^-function regularization, after explicitly 
computing the spectra. 

But this is not the whole story: As we are going to see, these global boundary conditions give rise to the presence 
of zero modes, which must be taken into account. 

III. ZERO MODES 

We will here show that the operator (-Dq)k has |fc + 1| zero modes. 
From (0), we get 

D„ i/; = 0^^e-"'^='^('') 7/' = 0, (16) 

or, equivalently 



e-*«(a,. - iSe) \ / e-^-^M \ ( Lp(r,6) 



0. (17) 



Now, we introduce the expansions 



n—-'OQ 
oo 



n— — OQ 



The solutions are thus given by 

Xn{r) = bn r-" e-"^('-), (18) 

where the coefficients a„ and 6„ are to be determined from the nornializability requirement at the origin and the 
boundary conditions at r = i?, Eq. (O). Thus 
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^{r,6) = e"^M^ an r" e^"^ 

n=0 

e-"*(^) ^ 6; r" e-'"^ (19) 



n=0 

Then, there are |fc + 1| zero modes which, once normahzed, are given by 

ga0(r) / 



V27r g„(i?;a) V 

-Q0(r) / Q 



for < n < fc, if fc > 0; 



Here, the normalization factors are 



, for < 71 < -A: - 2, if fc < -1. (20) 



Ja 

Pn{u;a) = / e-^^-^W r2"+i dr . (21) 



Notice that, for fc = — 1 (in particular, when $ = 0), there is no zero mode. 
So, the kernel of the orthogonal projector on Ker(_DQ,) is 



-fe-2 



g-Q[0(z) + 0(lu)] / Q Q 

P„(z,u;) - 2^ n ^^.T^y. , if fc < ~1. (22 



n=0 



2^p„(i?;a) VO (^*^) 



Since Pa is an orthogonal projector, = Pai we have 

9aPa (1-P„) =Pa S^^a- (23) 

Being (Z3q + Pa)K invertible, we can define 

Det{Da)K = Det{Da + Pa)^, (24) 

and write 

Det{D)n _ Det{D + Pi)n Det{t ^ + Po)k 



Det{i^)K=o Det{i^ + Po)K Det{i ^)k=o 



As mentioned above, we can compute 

d 



(25) 



„ ,lnDet{Da + Pa)K,]^Tr[i4 + dMGa], (26) 
oa 

after determining the Green's function G{x, y) satisfying 

(Da + Pa) G{x,y)^5{x,y), 

{V> V<)Gix,y)\r^n^O, (27) 

and applying a regularization prescription to the trace in ( p^ ) . Then, by integrating in a from to 1 we will get the 
first quotient in the r.h.s. of (p5|). 

The computation of this Green's function will be the subject of the next section. 
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IV. CALCULATION OF THE GREEN'S FUNCTION 

Notice that, since Da is self adjoint, we can write 

G(x,y) = (1-P„) (1-P„)+ Pa, (28) 

where Q{x,y) is the kernel of the right-inverse of Da on the orthogonal complement of Ker{Da)- From the Eq. (|^), 
it is easy to see that 

g{x,y) = 6"-^^*^ ^0(^^,2;) e"T='^('-'), (29) 
where Go{x, y) can be obtained by solving the differential equation 



for 



satisfying 



Pa e"''^'^^'^) Vx) " ^' 
{V> P< )e«^=^W (^^) =0. 



(31) 



Now, by expanding 



^(r,0)=^<^„(Oe^"', Xir^e) ^Y.^n{r) 
f{r, 9)^Y1 0) =J2 9n{r) e™^ (32) 



we get 



n>fc+l ''^ n<k •'^ 

f Tdr' rde' e™(«-«')(^)"e-'^g(-',^')-f T T d^' J] e™(''-'^') (^)" e-«'g(r', 0')- (33) 



n>fc+l n<fc 

Analogously 



/•-K /•27T / ' \ " ' f-^''^ / ' \ " 

X(r,0) = f / dr'/ d^?' ^ e™('^-'^')(-) e^^' f{r' ,9') - ^ dr' dQ' e'^^'^''^ [-] e^'' f{r' ,6'). 

^■^Jr Jo „<fc+i \rj ZTT Jo Jo 

n>k+2 ^ ' 

(34) 



Therefore, we have 



1 / 1 



^o(x,.) = -|^^^^..i "0 J' ^''^ 

and 



1 / - 7^7"^' (f ) 



Finally, from Eq. (28), after a straightforward but tedious computation, we obtain, for fc > 0, 
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/ fc 



oa[0(a;)-0(y)] 



k 

1 I J. V /'2L'\" 
X-Y Y ^\ y) g„(„., 
n— 



g-a[0(a;)-0(j/)] 



fe 

E 

n=0 



X'-Y* X* l^^KX'l q^iR.a) 



(37) 



and, for /c < — 1, 



-a[0(a;)-0(y)] 



X 





-fe-2 

*-Y* ^ Y' ^ \Y' 
n=0 



„a[ij>(x)-(t>{y)] 



X 



1 L (^\" P^i^:" 

-Y X 2^ \X) p^^a:a 

n—0 



^g-a[0(2:)+0to)] y- jX'Y)" 



(38) 



V. EVALUATION OF THE DETERMINANT 



Now we have all the necessary elements to compute the first quotient in the r.h.s. of Eq. (p5[). In fact, from Eqs. 
(H) and (H) 



rr[(9„P„)G„] =0. 



(39) 



So, Eq. (E6[) becomes 



dc 



[lnDet{D„ + P^)^]^Tr[4Go 



(40) 



As usual, the kernel of the operator inside the trace is singular at the diagonal, so we must introduce a regularization. 
We will employ a point-splitting one where, following Schwinger we will introduce a phase factor in order to 
preserve gauge invariance. We thus get, for > 0, from ( p7| ) 

4(a;)G„(a:,a; + e)e"^-^(^) 



Tr [4Ga] = / d'x tr 



2ti 



r<R 



k 

-T 



Qn{x;a) 

€r - iee ■ X* qn{R;a) 



0{e) 



-ie 



E 



0(e) 



+0 "'O 



dr r(\) (r) 



a 



r<R 



d^x 0'' -2{k + l) cj){R) + ^ — In [g„(i?; a)] 



n=0 



(41) 



Finally, performing a "symmetric" e — > limit, which drops out the first term, and integrating in a from to 1, 
we get 



In 



Det{D + Pi)K 
[Det{i ^ + Po)k 



I d^x 0'^ - 2 (fc + 1) 0(i?) + ^ In 



27r 



r<R 



2{n+l) 



qn{R;i) 

i?2(«+l) 



(42) 



When there are no zero modes (fc + 1 = 0) only the first term in the r.h.s. survives. For k < —1, from (p8|), one 
gets a similar expression where the sum runs from to — fc — 2 = |fc + 1| — 1, and the norms, defined in Eq. ([2l]), 
Qn, are replaced by Pn- This result can also be obtained by performing the change cj) ^ —(f) and fc + 1 — (fc + 1). 
Moreover, as expected from (^, this expression is invariant under the shift (j>{r) — > (/)(r) + e, for any constant e. 

In the following, we will obtain the second quotient of determinants in Eq. (|2^) by computing explicitly the spectra 
of the free Dirac operators and using a ^-function regularization. The equation to be solved is 
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or, explicitly 



Now, by expanding 



i ^ij) = Xip, 
e-''{dr-^de)\ (fir.e) 

oo 

n— — OO 

OQ 



Xir,0) 



(43) 
(44) 



(45) 



71 — — OQ 



we get 



n A 

« <^U^) - i -Vn{z) = Xn+l{z) 



^ Xn+l{z) + * Xn+l{z) = (^ra^j 



(46) 



where z = |A|r. So, </Jn(z) satisfies 



¥':(^)+ -(^;(z) + (i-^)v'n(^) = o. 



Therefore, the general solution of Eq. (^4|) is 

oo 

7/;(r,i7)-^^„ 



n=0 



Jn(|A|r) 
-z^J„+i(|A|r) e^("+i)e 



(47) 



(48) 



Now, by imposing the global boundary conditions given in Eq. (|1 1|) , it follows that: If n > fc + 1 , the upper 
component must vanish at the boundary, which implies A = ±jn,i/R ( jn,i is the l-th zero of Jn(z)). Analogously, if 
n < fc, it is the lower component which must vanish, and so A = ±j„+i,;/-R- Therefore, the eigenvalues are 



Xr, 



±jn.i/R, for n = 0,±1,±2, ... and 1 = 1,2, 



(49) 



Notice that j-n,i = jn,i, and that, for n = k + 1 the eigenvalues appear twice, once for an eigenfunction with vanishing 
upper component at the boundary, and once for another one with vanishing lower component. 
From these eigenvalues we can construct the ^-function 



n— oo oo 



[_n— — oo l — l 



1=1 



J|fc+iM 
R 



(50) 



The first term, |fc + 1|, is the multiplicity of the 0-eigenvalue of {i P)k- It is also interesting to note that the double 
sum in the r.h.s. corresponds to the ^-function of the Laplacian on a disk with Dirichlet boundary conditions, thus 
being analytic at s = . The second sum, to be explicitly computed below, is regular at s = 0. Then C{i^+Po)K{^) 
is regular at the origin. As far us we know, the (^-regularity at the origin for non local boundary conditions has not 
been established in general 
From (|^), we can write 



In 



Detji ^ + Po)k 
Det{i 



lK=0 



d_ 

ds 



(1 + e— ) T ^( - J2 Oo^')" 



. 1=1 



1=1 



s=0 



/|Wi|(0) - /o(0) + (Ini?- y)[/|..+i|(0) - /o(0)] 



(51) 
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where the function / is defined as 



1=1 



Taking into account the asymptotic expansion for the zeros of Bessel functions ||T^ 



( J,,)-' - (In)-' + s 



2v-l 



T:{lTr) 



we can write 

oo 



(=1 



(52) 



(53) 



(54) 



where C,{s) is Riemann's C-fmiction. Notice that the first series converges for Re s > —1; thus, it can be evaluated at 
s = 0, to obtain 



MO) 



V 1 

2^4' 



and 



/:(0) ^ ln2 + f ) (luTT - 7) - 5] hi 

1=1 



-— e V 4 1 y 

ITT 



2 V 4 

where 7 is Euler's constant. By replacing Eqs. (fc5|) and (|56|) into Eq. (fell), we get 



(55) 



(56) 



In 



Det{i ^ + Po)k 



L Det{i ^)k=o 



+ -7-ln(^)] + 25:in 



1=1 



Jo,i 



(57) 



Finally, adding ( }42D and (|57|), and taking into account that we have used a gauge invariant procedure, we obtain 
for the quotient in the l.h.s. of (P5|) 



In 



DetiD + Pi)K 



El- 



2{n + 1) 



Det{i ^)k=o 

9n(i?;l) 



__L /■ rf2^ A^(5^^_M!i) A,-2 (fc + 1) 0(i?) 



i?2(n+l) 



|fc + l|[^-7-ln(-)]+2Eln 

1=1 



(58) 



The first term is the integral on the disk of the same expression appearing in the well-known result for the boundaryless 
case po| . Note that this final result is also invariant under the transformation: — s- + e for any constant e. 

Now, a comment is in order concerning global axial transformations and their relationship to zero modes. Under 
such a transformation, with constant e. 



e-'^^t ^)n=^e-<^'' = {^ ^)^=o 

is invariant. Moreover 

e-75^p + Pi)Ke-^=^ ^{D + e-T^^Pie-'^^^)K, 
while its inverse has the structure 

G(')(a:,y) = (1 - Pi) G^x^y) (1 - Pi) + g-r^^ p^e^^' 
Therefore, since 75 leaves Ker D invariant, 



(59) 
(60) 
(61) 
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d_^^ I Pet (e-T^^p + PQ^e^T^^) 
de [ Dei (e-T5<^(i ^)K=oe-^'3') 



= -Tr 



{75, Pi} e-^='G(^)l = -2Tr [75F1] = -2(iV+ - iV_), 



(62) 



where -^+(-) is the number of positive (negative) chirahty zero modes. 

It can be verified that our strategy leads to the correct resuh for the index of D. By foUowing the same procedure 



that lead to Eq, 
instead of Eq. (|2S 




58|), we can compute the quotient of determinants in the l.h.s of (|6 
the only difference appears in the first term of the r.h.s. of ([50|), 



^). In fact, using Eq. (|6l| ) 
where a factor e^^'^'^ arises. 



Thus, after performing the e-derivative 



7V_ 



1, 



which agrees with our previous result for the number of zero modes. 

For the sake of completeness, in the next section we will show how the presence of the ct's in 
the index theorem. 



(63) 



and (O) affects 



VI. RELATION TO THE INDEX THEOREM 



Here we will reobtain the index by following the Atiyah-Patodi-Singer approach |13| 
The Dirac operator in polar coordinates can be written as [15 1 











(64) 



subject to the boundary conditions in Eq. (^ij). It is to be stressed that the presence of the factors e^**/^ coming 
from the cr's will produce a departure from the index formula in [p^ . 

To be be consistent with (qTI) we call V the differential operator {~dr + B) acting on functions 



/(r,0)=^/„(r)e^("+^) 



satisfying 



V> 



-iB/2 



Its adjoint, X>+, is the differential operator {dr + B) defined on the space of functions 

n 

5(r,0) = ^g„(r)e*("-^)« 



satisfying 



Then D+D 



= 0. 



i9/2 25+25 gi9/2_y^^ -^^j^gj-g 25+25 = —Q"^ -|- 52 acting on functions satisfying 

V< [e'''/'^Vf{R,e)\ = 0. 



or, equivalently, 



fn{R) = 



for 71 > fc + 1 



-f:^{R) + ''^^'l "" UR) =0 for n < fc. 



(65) 



(66) 



(67) 



(68) 



(69) 



(70) 
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Similarly, D'D'^ = —d^ + acting on functions g{r, 0) satisfying 

5„(i?) = forn<fc + l 

(71) 

g'niR) + '-^-^^^^—^gniR) = o for n > fc + 2. 

Now, we have for the heat kernels 

Tr {e-*^+^} - Tr {g-*^^^ } = Tr {e-*^+^} - Tr {e-*^^+ } . (72) 

When written in this fashion, the boundary contribution to the r.h.s. in (fz^), K{t), can be easily constructed from 
Eqs. (2.16) and (2.17) in |l3), thus getting 

m = I eric (-{k +l-K)V-t)-Y: "'^"'"/"''^ erfc (| n + i - ) . (73) 



Taking into account that 



we obtain 



where 



In particular 



for k + ^ < K<k + 1 

lim K{t) = m{K) i for K = fc + i , (74) 

1 for k < K < k + ^ 



dt t^-' {Kit) - m(K)} = ^^l^^^^ ViB+^)i^s), (75) 
r/(s+^)(s) = ^ sig(n+i-«;) |r^+i-Kr^ (76) 



2{K-k-l) for fc+i<K<fc + l 
ViB+^)iO)=\ for K = . (77) 

2{k - k) for fc < K < /c + i 

Finally, taking into account that, ior k < K < k + 1, 

miK) - i ry(B+^)(0) = /c + l- K= i[l- /i(S) - ry(B)(0)] , (78) 
where h(B) = dim ker(i?), and from the asymptotic expansion of the heat kernels in (|7^), we have 

in agreement with (|63|). The first term in the intermediate expression is the well known contribution from the bulk 
po| . The second one is the boundary contribution of APS, shifted by 1/2. This correction, due to the presence of the 
factor (T in (|^), has already been obtained in with slightly different spectral boundary conditions. 



VII. CONCLUSIONS 



In this paper we have achieved the complete evaluation of the determinant of the Dirac operator on a disk, in the 
presence of an axially symmetric flux, under global boundary conditions of the type intoduced by Atiyah, Patodi and 
Singer iQ. To this end, we have proceeded in two steps: In the first place, we have grown the gauge field while 
keeping the boundary condition fixed. This calculation was possible thanks to the exact knowledge of the zero modes 
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and the Green's function (in the complement of the null space.) Here, a gauge invariant point splitting regularization 
was employed. In references ^0 11 1 we developed, for the case of local boundary conditions, a regularization scheme 
based on Seeley's complex powers Its application to the present problem leads to the same result as the point 
splitting we used here, even though its relation to the C-function is not guaranteed. In fact, as far as we know, the 
construction of complex powers for elliptic boundary problems with global boundary conditions is still under study 
@- 

In the second step, we have explicitly obtained the eigenvalues of (i ^+Po)k- We have shown that the corresponding 
C-function is regular at the origin and we have evaluated the quotient of the free Dirac operators for two different 
global boundary conditions. 

We have verified that our complete result is in agreement with the APS index theorem, adapted to our example. 
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